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Abstract. We review the current status of lattice QCD studies of the nucleon system. In particular, we focus on the 
determination of the shape of the nucleon by probing its wave function as well as by evaluating the N to A transition form 
factors. 



I. INTRODUCTION 

Lattice QCD provides an indispensable method in our efforts to solve Quantum Chromodynamics (QCD). It is the 
only approach, using directly the QCD Lagrangian, available up to now in the energy regime between high energy 
processes, where perturbation theory is applicable, and very low energies, where chiral perturbation theory can be 
used. Lattice QCD is a discretized version of QCD formulated in terms of path integrals on a space-time lattice 1 1 ] 
with only parameters the bare quark masses and the coupling constant, like the continuum theory. One recovers 
continuum physics by extrapolating results obtained at finite lattice spacing a to a = 0. In order to perform the 
continuum extrapolation a separate calculation at several values of a is required. Because this can require large 
computer resources, in most cases, the strategy is to work close to the continuum limit, either by choosing a sufficiently 
small or by improving the discretization scheme to eliminate order-a terms. Numerical evaluation of observables in 
lattice QCD necessarily requires that the size of the lattice is finite. To keep the size of the box large enough to fit 
the hadrons well inside one must increase the number of sites as one decreases a. Therefore calculations at decreasing 
values of a require increasingly larger computer resources. To keep finite volume effects small one must have a box 
that is much larger than the Compton wavelength of the pion. Usually we take Lm K ~ 5 where L is the spatial length 
of the box and m n the pion mass. In an actual calculation the bare quark masses for the u and d quarks are taken 
larger than in the real world. Given that computational costs increase like m~ 9 , the use of larger quark masses enables 
inversion of the fermionic matrix, needed for the calculation of hadronic matrix elements, with currently available 
resources. To extract physical quantities from lattice calculations typically done with pion masses above 350 MeV, 
one needs to extrapolate to the physical quark masses. It is only very recently that we can reliably reach pion masses, 
below 350 MeV \2i\M- This is a very important step forward in our effort to eliminate one source of systematic error 
associated with the extrapolation to the light quark masses. Typical state-of-the art lattices use a ~ 0. 1 fm and L ~ 3 fm 
for pions of mass down to about 350 MeV. 

Discretization of space-time introduces an ultra-violet cut-off limiting the highest momentum to 2n/a. Therefore, 
the lattice provides a regularization of the ultra-violet divergences of the theory making it well defined from the start. 
Working in a finite box allows only discrete values of momentum in units of 2j[/L. The bare coupling constant and 
quark masses are tuned as a changes to leave physical quantities unchanged. Besides using a finite lattice spacing, 
a finite volume and heavier u- and d- quarks, an additional step, that enables us to numerically evaluate the path 
integrals needed, is rotation to imaginary time, t — > —it, resulting in replacing exp(/5) by exp(— Se) where S (Se) 
is the classical QCD action in Minkowski (Euclidean) space. Having a real action one can apply stochastic methods 
commonly used in statistical mechanics to evaluate the path integrals. Whereas finite a, volume and quark masses 
are amenable to systematic improvements, rotation to Euclidean space selects a set of observables that can be studied 
within this framework. These are observables that can be determined from the properties of the discrete lower lying 
states. For this set of observables, lattice QCD produces the exact answer provided the extrapolations to the continuum 
and infinite volume limits as well as to the physical quark masses are under control. 

In a typical lattice calculation one starts by choosing the bare coupling constant g, which fixes the lattice spacing, 
and the bare masses for the u-, d- and s-quarks. One then computes a physical quantity such as the mass of the pion 



and the nucleon in lattice units as a function of the quark mass. The pion mass is used to fix the u- and d- quark 
masses (assumed degenerate) and the mass of the kaon or <j> to fix the strange quark mass whereas the lattice spacing 
is determined by extrapolating the results, for instance, for the nucleon mass to the physical pion mass. Any other 
physical quantity in the light quark sector then follows. If instead of the nucleon mass one chooses another physical 
quantity such as the pion decay constant, fa, to set a the resulting value should be the same if finite lattice volume and 
non-zero a effects are under control. 

Like in the continuum, the lattice QCD action, Se = S g [U] +Sf[U, yf, y/], has a purely gluonic part S g written in 
terms of the gauge link U^(n) — exp(/a^A| i (n)), which connects site n with n + 1 in the /x-direction and a fermionic 
part Sf, which contains the kinetic energy of the quarks and the interaction terms. (n) denotes the gluon field. Gauge 
symmetry is exactly preserved by Se- As in any renormalizable theory, the fermionic action is bilinear in the fermion 
fields and can be written in terms of the fermionic matrix D as Sf — L„._/ W{ n )DnjV(j)- The exact form of D depends 
on the discretization scheme used for the fermions. The simplest is due to Wilson 1 1]. It has been widely studied but 
has the disadvantage of breaking explicitly chiral symmetry. Recent theoretical developments made it possible to have 
chiral fermions on the lattice. There are two equivalent formulations known as domain wall fermions ||J,|5|] and overlap 
fermions |6, 7]. They both require larger computer resources than Wilson fermions. 

The vacuum expectation value of any gauge invariant operator O can be computed by evaluating the path integral 

1 1 Jd[U}d[f]d[\ir]e- s siu}-SF[u,vM ' <• ; 

Integrating over the fermionic degrees of freedom we obtain 

[d\U]det(D\U])0\U,D- l \U]]e- s *M r 
<Q.\0\Q.>=— — z — Z = d[U] 6et(D[U})e-^ [u] , (2) 

where a factor DJ*[U] substitutes each appearance of —\j/ n \f/j in O. One can now perform the path integrals nu- 
merically by stochastically generating a representative ensemble of gauge fields U according to the probability 
exp { -S g [U] + In (det(D [£/]))} /Z and then compute 

1 N 

<£l\0\Q. >= lim - V 0[U k ,D- l [U k }] , (3) 

which involves the evaluation of the inverse of the fermionic matrix. For a typical lattice of size 24 3 x 48 the dimension 
of the complex matrix D is 8 million by 8 million. Therefore the time consuming part of a lattice calculation is the 
generation of an ensemble of gauge configurations and the computation of the inverse of the fermionic matrix D which 
yields the quark propagator. In many applications only a column of D 1-1 is required. 

The first lattice calculations were performed in the quenched approximation, which neglects pair creation by setting 
det(D) = 1 in Eq. (|2jl. This facilitates the generation of gauge links since one is left with the local action S y [I/]. 
In the quenched approximation the program outline above, i.e. taking the continuum and infinite volume limits and 
extrapolating to the physical quark masses, has been carried out for the spectrum of the low lying hadrons 1 8]. During 
the past five years theoretical progress in combination with terascale computers have made unquenched calculations 
with light pions and large enough volumes feasible using a number of different discretization schemes. 

In order to study the role of the pion cloud, which is expected to provide an important ingredient in the description 
of the properties of the nucleon system, one must generate dynamical gauge configurations with light quarks using 
large volumes. In this work the light quark regime is studied in two ways: 

1 . We use configurations with the lightest available dynamical Wilson fermions. The unquenched configurations are 
simulated with two degenerate flavors of Wilson fermions flolfTlll . 

2. We use MILC configurations generated with two degenerate light and one strange staggered quarks using the 
Asqtad improved action |2]. For the valence quarks we use domain wall fermions that preserve chiral symmetry 
on the lattice. This is therefore a hybrid calculation that uses different fermions for the sea and valence quarks, 
i.e the matrix D appearing in the determinant in Eq. @ is different from D^ 1 involved in the calculation of O in 
Eq. 0. Such hybrid calculations have been successful in recent evaluations of fundamental physical quantities 
such as gA 1 12] and the pseudoscalar decay constants fa and fa II 311 . 



Bearing in mind that both quenched and unquenched Wilson fermions have discretization errors of order a, and that 
both Asqtad and domain wall fermion actions have discretization errors of order a 2 and fermions preserving chirality, 



in contrast to Wilson fermions, agreement between the results within these two different lattice fermion formulations 
provides a non-trivial check of consistency of the lattice results. The hybrid calculation is computationally the most 
demanding since it requires propagators on a five-dimensional lattice. The bare quark mass for the domain wall 
fermions, the size of the fifth dimension and the renormalization factors Zy and for the four-dimensional vector and 
axial vector currents are taken from Ref. . The parameters of our calculations are given in Tabled 

TABLE 1. We give the number of configurations, the hopping parameter, K, which determines 
the bare quark mass via the relation 2am q = ( 1 / ic — \/k c ) for the case of Wilson fermions or the 
mass of the u and d quarks, m/, for the case of staggered quarks, the pion, nucleon and A mass in 
lattice units with their statistical errors determined from a jackknife analysis. For Wilson fermions 
a is set using the nucleon mass at the chiral limit whereas for staggered fermions we take the value 
extracted from the static qq force as determined in Ref. fl4ll . 



no. of configurations 


K (Wilson) or ami (staggered) 


am K 




a Ma 






Quenched 32 3 x 64 a~ l = 


2.14(6) GeV 




200 
200 
200 


0.1554 
0.1558 
0.1562 
Kb =0.1571 


0.263(2) 
0.229(2) 
0.192(2) 
0. 


0.592(5) 
0.556(6) 
0.518(6) 
0.439(4) 


0.687(7) 
0.666(8) 
0.646(9) 
0.598(6) 




Unquenched Wilson 24 3 x 40 a 


1 =2.56(10) GeV 




185 
157 


0.1575 
0.1580 


0.270(3) 
0.199(3) 


0.580(7) 
0.500(10) 


0.645(5) 
0.581(14) 




Unquenched Wilson 24 3 x 32 a 


1 =2.56(10) GeV 




200 


0.15825 
(Q- = 0.1585 


0.150(3) 
0. 


0.423(7) 
0.366(13) 


0.533(8) 
0.486(14) 


MILC 20 3 x64 «"' = 1.58 GeV 


150 
150 


0.03 
0.02 


0.373(3) 
0.306(3) 


0.886(7) 
0.800(10) 


1.057(14) 
0.992(16) 


MILC 28 3 x 64 a' 1 = 1.58 GeV 


118 




0.01 


0.230(3) 


0.751(7) 


0.988(26) 



The main goal of this program is to calculate within lattice QCD the fundamental physical quantities of the nucleon- 
A system. Providing the complete set of form factors and coupling constants constitutes a very important input for 
model builders and for fixing the parameters of chiral effective theories. In this presentation we discuss our results for 
the nucleon elastic form factors and the electromagnetic and axial N io Is. transition form factors as well as show first 
results on their wave functions. Other recent lattice studies on nucleon properties can be found in Ref. 1 15]. 



II. LATTICE TECHNIQUES 

The vacuum expectation value of gauge invariant operators is computed by numerical evaluation of appropriately 
defined path integrals. Let us first consider the evaluation of hadron masses. The vacuum expectation value of the time 
ordered product G h {t,q) —< £2|£ x exp (z'q.x) ffhfat)/} (0)\Q. > can be evaluated using Eqs. (0 and which require 
inverting D for each non-degenerate quark flavor once per gauge configuration. The interpolating fields, Jf,(x) are 
operators in the Heisenberg representation that create a trial state with the quantum numbers, h, of the hadron that we 
want to study. For example in the nucleon case an appropriate interpolating field is Jn{x) = £ a bc [u a T (x)CYsdb(x)] u c (x) 
where C is the charge conjugation operator and Latin indices denote color quantum numbers. The large time behavior 
of the correlator G (t ,0) yields the mass: 

n,p x n 

G h (t,o) '=4? i <a\j h \h ,o> |V"V (4) 



where E n (q) = ^rn^ + q 2 , \ho > is the lowest eigenstate of QCD with quantum numbers h with mass m/ 1Q and we have 
taken q = 0. For large time separation t between the source and the sink the unknown overlap factor | < \ho, > | 2 



and exponential dependence cancel in the ratio m e ff(f) = — log \G h (t)/G h {t — 1)], which therefore becomes time 
independent and can be fitted to a constant to extract the mass, m/, , of the lowest state. 

Three-point functions of the form G M ' ! (?,?i;q) =< £2|L x , y e i(l - x ff h (y,t)A(x,ti)J^(0)\Q. > are required for the 

evaluation of form factors. Inserting a complete set of hadronic states between operators in G llAh (t ,fi;q) as we did 
above for the extraction of hadron masses from two-point functions, we obtain 

G hAll {t,t x ;q) = £ £ e iqx < «|/»|«,p' > e -W)yp'.y e £»(p'h e -ip'.x < n rf\A\k,p > e -E n {v)h e 0<x < jt,p|JJ|Q > 

A ; ;i,p',px.y 

A £ (r-,)>>UA £(l »i < a]fhlho ><c ~ hoM ft h p >< h ,0\Mho,P > e- m ^e-% iph (5) 

where AE and AE are the energy differences between the two lowest hadronic states with quantum numbers h and h 
and p = — q. The exponential time dependence and unknown overlaps can be canceled by dividing with appropriate 

combinations of two-point functions G h (t,q). For example in the ratio R = G hAjl (t , t \ ; q)/ \J G h {It - 2t\ , 0)G S (2t x , q) 

the overlap factors and exponentials cancel in the large time limit when the ground states /io and Iiq dominate, yielding 
the matrix element < /io,0|A|/io,p >■ Since both two- and three- point functions decay exponentially it is crucial in 
constructing the ratio to choose combinations of two-point functions that involve the shortest possible time separations 
and to use techniques that isolate the lowest hadronic states \Iiq > and ho > at short time intervals. The former can 
be don e by choosing better but more complicated ratios than the one given here, to be discussed in Sections IV and 
V I16lll7il . whereas the latter by constructing better interpolating fields. Smearing techniques are routinely used for 
achieving ground state dominance before the signal from the time correlators is lost in the noisy large time limit. We 
use gauge invariant Wuppertal smearing to replace local by smeared quark operators at the source and the sink 111 8ll 
and, when needed, we apply hypercubic averaging 1 19] of the gauge links that enter the construction of the Wuppertal 
smearing function. As can be seen from Eq. l|5}, to compute three-point functions two sums over the spatial volume 
are needed. Performing the Wick contractions on the quark level for the three-point function one finds expressions 
that involve the full inverse of the fermionic matrix D (all-to-all propagator). This is to be contrasted with two-point 
functions for which only one spatial sum enters and therefore only one column of D is needed. The technique to 
automatically do one of the spatial sums by using an appropriately defined input vector that uses one column of the 
inverse of D when inverting, is known as sequential inversion. One can choose which of the two sums to do first. In 
this work we do the sum over the sink y. This requires that the quantum numbers of the sink are fixed but allows 
any operator A with any momentum q to be inserted. This means that with one sequential inversion one can extract 
the matrix element < liQ,0\A\ho,p > for given hadronic states but different operators and momentum q. We therefore 
measure the matrix element for all lattice momentum vectors that result in a given momentum transfer squared q 2 , 
thereby obtaining many statistically independent evaluations of the same form factors reducing statistical noise. 



III. PROBING HADRON WAVE FUNCTIONS 

Our main motivation for studying density-density correlators is that they reduce, in the non-relativistic limit, to the 
wave function squared yielding detailed, gauge-invariant information on the internal structure of hadrons 120112111 . The 
shape of hadrons, which is the topic of this workshop, is one such important quantity that can be directly studied. 
There are indications from experimental measurements of the quadrupole strength in the y*N — > A transition 112 2ll 
that the nucleon and/or the A are deformed. Density-density correlators, shown schematically in Fig.^ are four-point 
functions since they involve the insertion of two density operators. Four-point functions are technically harder to 
compute than three-point functions discussed in the previous Section. In particular, they require computation of all the 
spatial columns of the inverse of the fermionic matrix requiring L 3 inversions. A straight forward computation of such 
an inverse is therefore prohibitively expensive. In our first study I23B . which we consider as a feasibility study, the 
density-density correlators were evaluated without explicit projection to zero momentum hadronic states. This means 
that higher momentum states are suppressed only by the Euclidean time evolution. 

We have analyzed 220 quenched configurations at /3 =6.0 for a lattice of size 16 3 x 32 obtained from the NERSC 
archive, using the Wilson Dirac operator with hopping parameter K = 0.15,0.153,0.154 and 0.155. Using the relation 
2am q — (1 / K" — 1/jc c ), with the critical value K c = 0.1571, one can obtain the naive quark mass m q . A physical 
dimensionless quantity that is sensitive to the value of the bare quark mass is the ratio of the pion mass to the rho mass, 
which at these values of K, is 0.88, 0.84, 0.78 and 0.70 respectively. We fix the source and the sink for maximum time 




\ (r'+r,|) / 



D(r,t u t 2 ) = J d\' {h\p u {r' + r,t 2 )p d {r',H)\h) 

pt{x,t) = :/(r,f)7o/(r,f): (6) 

( x o>^) where for the density operator we use the normal order 
product so that disconnected graphs are excluded. 



FIGURE 1. The density-density correlator, D^c,t\,tj), for a meson. The time separations t\, Pi, t — t\ and t — t 2 are taken large 
enough to isolate the mesonic ground state. 





2.5 



7.5 



10 



FIGURE 2. A three-dimensional contour plot of the correlator is shown in red. Left: for the rho state with spin projection (cigar 
shape); Right: for the A + state with +3/2 (slightly oblate) spin projection for two dynamical quarks at K = 0.156. Values of the 
correlator (0.5 for the rho, 0.8 for the A + ) were chosen to show large distances but avoid finite-size effects. We have included for 
comparison the contour of a sphere (green). 



separation, which for this lattice is tja = 16 given the antiperiodic boundary conditions in the time direction. Both 
density insertions are taken at the same time slice at the middle of the time interval. To investigate the importance 
of dynamical quarks, we use SESAM configurations |24] generated with two dynamical degenerate quark species 
at j3 = 5.6 on a lattice of size 16 3 x 32 at K = 0.156 and 0.157. The ratio of the pion mass to rho mass is 0.83 at 
K = 0.156 and 0.76 at K = 0.157. These values are close to the quenched mass ratios measured at K = 0.153 and 
K = 0.154 respectively, allowing us to make pairwise quenched-unquenched comparisons. The general conclusion of 
this comparison is that, at these heavy quark masses, unquenching effects are small. We show in Fig.|5]contour plots 
of the density-density correlator for the rho and the A for the case of dynamical Wilson fermions at the heavy quark 
mass. The elongation in the rho is clearly visible whereas the A + appears to be squeezed. 

We have recently developed techniques for the evaluation of density-density correlators with explicit projection to 
zero momentum hadronic states l25ll and we are in the process of analyzing unquenched SESAM configurations at 
lighter quark masses as listed in Tabled 



IV. NUCLEON FORM FACTORS 



Recent results from polarization experiments 

Elm, 

have shown a qualitative different behavior than the tradi- 
tional Rosenbluth separation for the ratio of the proton electric to magnetic form factor, pL p G p E /G p M . An accurate 
determination of these form factors in lattice QCD can provide an important theoretical input for understanding the 
<7 2 -dependence of this ratio. We present here results obtained using Wilson fermions on the quenched and unquenched 
lattices given in Tabled The smallest g 2 -values are accessible on the quenched lattice, which has the largest spatial 
extent, since the smallest available non-zero momentum on a finite lattice is 2n/L giving — q ~ 0.17 GeV 2 . Although 
large momentum transfers are in principle available on typical lattices, the Fourier transform of two- and three-point 
functions becomes noise-dominated for momentum transfers beyond about 2 GeV 2 , limiting the range of high q 2 
values that can be extracted accurately. 

The standard decomposition of the nucleon electromagnetic matrix element for real or virtual photons is given by 



(N{p',s')\j^\N(p,s)) = 



M 2 N 



En(j?') e n (p) 



1/2 



u(p',s') 



2M N yH ' 



u{p,s) , 



(7) 



where p(s) and p'(s') denote initial and final momenta (spins), Mpj is the nucleon mass, F\ (0) = 1 for the proton since 
we have a conserved current and F 2 (0) measures the anomalous magnetic moment. They are connected to the electric, 
Ge, and magnetic, Gm, Sachs form factors by the relations 



G E (q 2 )=F l (q 2 )- 



(2M N ) 



:F 2 (q 2 ) G M (q 2 )=F l {q 2 )+F 2 (q 2 ) . 



(8) 



The electromagnetic matrix element is extracted from the three-point function G N ^ N following the procedure outline 
in the Section II. We use the lattice conserved electromagnetic current, 



E it 

f=u,d 



yj/f(x + fi)(l + J^ul^w 3 \x) - W f (x){l - y^U^yfix + n) 



(9) 



symmetrized on site x by taking jf (x) 
behavior of the improved ratio 



[j^(x) + j^(x — jl)] /2. We look for a plateau in the large Euclidean time 



^(f,f!;q;r;^) = 



<G^( f , fl ;q;r)> 



(G N (t,o-r 4 )) 

A£(r-fi)»l,A£;i»l 



(G N (t 2 -t u p;T 4 )) (G N (t h 0;T 4 )) (G N (t,0;T 4 )) 1 1/2 



(G N (t 2 -t h 0-r 4 )) {G N (t uP ;T 4 )) (G N (t,p-r 4 )) 



n( q; r ;M ) 



where we show explicitly the dependence on the projection matrices r, = | 



(10) 

j = 1,2,3 and T4 = 



for the Dirac indices. AE is the energy difference between the nucleon and its excited P\\ -state and 

2 



jj 


/ 


q = p for a nucleon at rest in the final state. Since we want to study the g z -dependence of the form factors we 
evaluate the three point function with sequential inversion through the sink. We fix the source-sink time separation 
t j a = 11(12) for the quenched (unquenched) Wilson lattices and search for a plateau of R(t ,fi;q as a function 

of the time slice, t\, at which couples to a quark. Q 2 = —q 2 denotes the Euclidean momentum transfer squared. 

The two Sachs form factors are extracted from the ratio defined in Eq. (II Oi by choosing appropriate combinations 
of the direction of the electromagnetic current and projection matrices F. Provided the Euclidean time separations t \ 
and t — t\ are large enough to filter the nucleon ground state, the ratio R(t,t\\q becomes time independent and 

the two form factors can be extracted from the Euclidean space expressions, 



£n(q;ry,M = 

k=i 



K 

2M N 



(P2-P3)(>lj+{P1-Pl)(>2,i + (PI -Pl)&,j \G M {Q 

n(q;r 4 ; ii = i)=K^- G E (Q 2 ) 



~2M N 



U( q ;T 4 ^=4)=K^±^G E (Q 2 ) 



2M N 



(11) 



(12) 
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FIGURE 3. The isovector form factors, Ge, (right) and Gm, (left) as a function of Q 2 . By Nf = we denote quenched 
lattice results at K = 0.1554 (crosses), at K = 0.1558 (open circles) and at K = 0.1562 (asterisks). Results using two degenerate 
flavors of dynamical Wilson fermions are denoted by Nf = 2 at K — 0.1575 (open triangles), K = 0.1580 (filled circles) and at 
K — 0.15825 (open squares). The filled triangles show experimental results for the isovector form factors extracted by interpolating 
the experimental data for the proton and neutron form factors (the details of the interpolations)are given in Ref. 1 17|l). 



where K = 



2M-, 



e n {e n +m n ) * s a f actor due to the normalization of the lattice states. We note that the expression for 
Gm is obtained by using an optimal linear combination for the nucleon sink that provides the maximal set of lattice 
measurements from which Gm can be extracted requiring only one sequential inversion. Eqs. dl 2i yield Ge with an 
additional sequential inversion. 

The Y*N — ► jV transition, in addition to an isovector part, contains isoscalar photon contributions. This means 
that disconnected loop diagrams also contribute. These are generally difficult to evaluate accurately since the all- 
to-all quark propagator is required. In order to avoid disconnected diagrams, we calculate the isovector form factors. 
Assuming St/ (2) isospin symmetry, it follows that 



; p \( 2 -uy»u - l -dYd) \p) - { n^-ufu - l -dfd) \n) = { p | (ufu - dfd) \p). 



(13) 



One can therefore calculate directly the three-point functions related to the right hand side of the above relation, from 
which the isovector nucleon form factors 



G E (q 2 ) = G P E (q 2 ) - G" E (q 2 ), G M (q 2 ) = G P M {q 2 ) - G M (q 2 ) 



(14) 



can be extracted using Eqs. (^J and d!2l > by only evaluating connected diagrams. 

Besides using an optimal nucleon source, the other important ingredient in the extraction of the form factors, is 
to take into account simultaneously in our analysis all the lattice momentum vectors that contribute to a given Q 2 . 
This is done by solving the overcomplete set of equations P(q;/i) =A(q;/i) -F(Q 2 ) where P(q;jj.) are the lattice 

/ q e 

measurements of the ratio given in Eq. (II Oi having statistical errors Wf, and using the different sink types, F = I ^ 

and A is an M x 2 matrix which depends on kinematical factors with M being the number of current directions and 
momentum vectors contributing to a given Q 2 . We extract the form factors by minimizing 



(15) 



using the singular value decomposition of the matrix A. Given the fact that one can have a few hundred lattice 
momentum vectors contributing for a given Q 2 , the statistical precision is highly improved. There is an additional 



advantage arising from including momentum vectors q as well as — q in our analysis: The lattice conserved current 
given in Eq. (|9jl differs from the local electromagnetic current ijr(x)YpY(x) by terms of order a. However when we 
average over q and — q these order a terms vanish. 
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FIGURE 4. Left: The form factors F\ (upper) and Fi (lower) as a function of Q 2 at the chiral limit. The results extracted from 
experiment are shown by the filled triangles. Right: Chiral extrapolation of the magnetic moment (upper) and the r.m.s radii r\ 
(middle) and ri (lower). The solid line is the best fit to the effective chiral theory results. The dashed lines show the maximal allowed 
error band using the errors on the fitted parameters. In all graphs quenched results are shown by the crosses and unquenched results 
by the filled circles. Quenched results from Ref. |28] are shown by the asterisks. 

In Fig. we show quenched and unquenched results for the isovector electric and magnetic form factors using 
Wilson fermions at the three values of the quark mass given in Table [2 Both quenched and unquenched results 
decrease with the quark mass yielding a larger slope at small Q 2 , which is the expected behavior. However both 
quenched and unquenched results for Ge deviate more from experiment than Gm for the same quark mass. The two 
main uncertainties regarding lattice results are finite a-effects and the fact that the u- and d- quark masses are larger 
than physical. In the quenched theory one can envisage repeating the calculation using a finer lattice in order to check 
whether finite a-effects can explain this behavior. Evaluation of these form factors closer to the chiral limit with Wilson 
fermions is problematic and one will have to use other types of discretization for the fermions to study the dependence 
on lighter quark masses. To be able to directly compare with experiment, using the results of this work, we need to 
extrapolate Ge and Gm to the chiral limit. The quark masses employed in this work correspond to pion masses in 
the range 560 to 410 MeV in the quenched theory and 690 to 380 MeV in the unquenched theory. Pion cloud effects 
are expected to be small in this range of pion masses and therefore we expect a linear dependence of the results 



on m n . Our lattice data at these quark masses confirm this expectation and therefore, to obtain results in the chiral 



limit, we extrapolate the form factors linearly in m 2 . We show the linearly extrapolated lattice results for F\ and F2 in 
Fig.@] We observe agreement between quenched and unquenched results at the chiral limit. In addition, there is good 
agreement between lattice results for F2 and experiment, with perhaps small deviations at small Q 2 . This is not the case 
for F\ where the experimentally determined isovector F\ decays faster as compared to the lattice results. In a recent 
calculation, the quark mass dependence of the isovector magnetic moment and radii was determined. This was done 



within a chiral effective theory with explicit nucleon and A degrees of freedom 1 28, 29] . Therefore we can extract the 
relevant low energy constants and counterterms that enter in the expressions of these quantities in the effective chiral 
theory from fits to our form factors and obtain results at the physical pion mass. Fitting to our lattice data we obtain 
the curve shown by the solid line in Fig.@] The details regarding the fits are given in Ref. fnll . The dashed lines give 
the maximal error band determined by varying the values of the fitted parameters by their errorbars. The extrapolated 
value of the magnetic moment at the physical pion mass is in agreement with experiment. The resulting fits for the 
radii are also shown in Fig. |4] The pion mass dependence of the Dirac radius is not well reproduced. Since this is 
related to the slope of F\ this is not surprising given that the lattice results have a different slope from the experimental 
results and hardly show any quark mass dependence. 



V. N TO A TRANSITION FORM FACTORS 

The evaluation of the N-A matrix element, within the fixed sink approach, requires a new set of sequential propagators 
since in the final state, instead of the nucleon, we have the A. However, once we produce the sequential propagators 
needed for the evaluation of the electromagnetic N-A matrix element, the axial one can be obtained with almost no 
additional computational cost using the same sequential propagators. The N to A transition involves no disconnected 
diagrams and therefore lattice results can be directly compared to experiment. 



V.l. Electromagnetic transition form factors 

To address the question of possible deformation in the nucleon system the experiment of choice is electroproduction 
of the A that measures the nucleon-A transition amplitudes. Non-zero quadrupole amplitudes are thought to be 
connected with a non-spherical nucleon or/and A 1 22] . We can establish direct contact with experiment by calculating 
the N to A transition form factors in lattice QCD. To obtain accurate results that can provide a meaningful comparison 
to experiment two novel aspects are implemented: 1) An optimal combination of three-point functions, which allows 
momentum transfers in a spatially symmetric manner obtained by an appropriate choice of the interpolating field for 
the A. This is similar to the construction of the optimal nucleon source for the calculation of the nucleon magnetic form 
factor, Gm, discussed in the previous Section but more involved 1 16]. 2) An overconstrained analysis using all lattice 
momentum vectors contributing to a given q 2 value in the extraction of the three transition form factors analogous to 
what was done for the elastic nucleon form factors. 

The matrix element for the y*N —> A transition with on-shell nucleon and A states and real or virtual photons has 
the form Q 

(A(pV) \f\N(p,s)) -'Vf ( ^^(p) )' 2 Ua(p',s')^u(p, S ) (16) 

where u a (p' ,s') is a spin-vector in the Rarita-Schwinger formalism. can be decomposed in terms of the Sachs 
form factors as 

^ = %1 (q 2 )K^ + &E2(q 2 )K^ + ^ci(q 2 )K^ , (17) 

where the magnetic dipole, S^n, tne electric quadrupole, @e2> an d the Coulomb quadrupole, ^cz< form factors de pen d 
on the momentum transfer q 2 — (p' — p) 2 . The kinematical functions K a ^ in Euclidean space are given in ref. [31]. 
The ratios Rem or EMR and R$m or CMR in the rest frame of the A are obtained from the Sachs form factors via 

p ^(<? 2 ) - |q| Sfc 2 (<? 2 ) n „ 

Rem — -jz — t-jt , Rsm — --z — tt, — 7-57 • U°) 



We look for a plateau in the large Euclidean time behavior of the optimized ratio 

(G A f N (t2A;q;r)) 



R a (t 2 ,h;q;r;n) 



■N»«u_ r/ G « (f2 _ fl5 p. r4)) ( G A(^0;r 4 )) (G$M;r 4 )) l 1/2 



(G£(? 2 ,0;r 4 )) 
To extract <$m\ we take the symmetric combination, 



(G$(fe-fi,0;r 4 )> (G N (t uP ;T 4 )) <G"(f 2 ,p;r 4 )) 



(19) 



Si(q;M) = L n CT (q ;r 4 ; M ) = iA £ e°^p^ Ml (Q 2 ) (20) 
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FIGURE 5. G*,(g 2 ) = l/y/l + Q 2 /(m N + m A ) 2 W M \{Q 2 ) as function of Q 2 . The notation for Wilson fermions is the same as 
that in Fig. [3] Results within the hybrid scheme are shown with the stars for m/ = 0.03, with the filled rhombi for m; = 0.02 and 
with the filled squares for mi = 0.01. Experimental results from Ref. 1331 are shown by the filled triangles. 
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FIGURE 6. Rem (left) and Rsm (right) as a function of Q 2 . Upper graphs show lattice results in the same notation as Fig. [3] The 
lower graphs show linear extrapolations to the chiral limit for quenched (crosses) and unquenched (filled circles) Wilson fermions. 
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dotted and dashed lines are the results from a dynamical model with bare and dressed vertices respectively II37I1 . 



computed for all spatial current directions so that lattice momentum vectors in all directions contribute. This com- 
bination, is built into the A interpolating field and requires only one inversion. Other combinations yield &E2 and 

In Fig.[5]we show our results for the magnetic dipole form factor at similar pion masses for quenched and dynamical 
Wilson fermions and in the hybrid approach using MILC configurations and domain wall fermions. Comparing results 
for quenched and unquenched Wilson fermions we see that unquenching effects are small in this range of quark masses. 
We also observe agreement between the results obtained in the two dynamical calculations. As we have already pointed 



out, agreement between the results using dynamical Wilson fermions and domain wall fermions is non-trivial since 
these discretization schemes have different lattice systematics. In all cases S^i decreases with the quark mass bringing 
lattice results closer to experiment. However, assuming a linear dependence in m\ to extrapolate lattice data obtained 
using Wilson fermions where the statistical errors are the smallest, we find results at the chiral limit, which are higher 
than experiment. Module finite a-effects, the conclusion is that to reconcile lattice results with experiment a stronger 
dependence on for small quark masses seems necessary. In the hybrid approach the errors are larger and a reduction 
in the errors by at least a factor of two is required in order to draw any concrete conclusions regarding whether they 
show a better agreement with experiment. 

Results for the EMR and CMR ratios are only shown for Wilson fermions in Fig.|6]since the errors on these ratios 
calculated in the hybrid scheme are too large prohibiting a meaningful comparison to the Wilson results. Although 
the results from dynamical Wilson fermions have larger errors than in the quenched case they tend to favor negative 
non-zero values for both ratios. However, given the size of the errors, an order of magnitude more configurations need 
to be analyzed in order to assess unquenching effects and draw a definite conclusion about the values of these ratios 
in the unquenched case. An analysis within a chiral effective theory leads to a non-trivial quark mass dependence for 
these ratios 1 40] that brings quenched lattice data in agreement with experiment at the lowest Q 2 -value. 



V.2. Axial from factors 



The N to A transition, besides being used to probe electromagnetic properties of the nucleon system, is also well 
suited for studying the weak structure functions. This is because the A(1232), as the dominant nucleon resonance, can 
be well identified and being a purely isovector spin-flip transition, provides selective information on hadron structure. 
A lattice calculation of the axial form factors is timely and important given that experiments at Jefferson Lab |41J] are 
underway to measure these form factors. In this work we evaluate the dominant contribution to the parity violating 
asymmetry, which is determined by the ratio C5 /C3 and is to be measured at Jefferson Lab. This ratio is the off- 
diagonal analog of the gA/gv ratio extracted from neutron j3 -decay and therefore tests low-energy consequences of 
chiral symmetry, such as the off-diagonal Goldberger-Treiman relation. In addition the ratio of axial form factors 
Cg /C5 provides a measure for the conservation of the axial current |42|]. 

The invariant N to A weak matrix element can be expressed in terms of four transition form factors as l43H44ll : 



<A(p',s')\Al\N(p,s)> = = 



M A M A 



\ 1/2 
] u 



3 \E A (p')E N (p) J 



C A s{q 2 )g X , ' W 



M 2 N 



V,*0 



u{p,s) 



~jj^~ 7 + ~mT P )\SX l igpv-gXpg ii v)q' J 



(21) 



where A 3 (x) = \^{x)y il ys^'^f{x) is the isovector part of the axial current (t 3 being the third Pauli matrix). Having 
evaluated the electromagnetic N to A transition form factors on the lattice this matrix element can be evaluated 
without requiring any further inversions since the optimized A sources are the same as those used in our study of 
the electromagnetic N to A transition and only the operator that couples to a quark line differs. We consider the same 

ratio as that given in Eq.Q9]but replace the three-point function (t ,h;q;,T) with (t ,fi;q;T). In the large 

Euclidean time limit this ratio yields the transition matrix element of Eq. (12 1 i . Using, for example, the A source, S\, 
for the evaluation of the three-point function we obtain for large time separations 1 1 and t — t\ 



Si (q, At -4) 



k=\ 



rA ,M A A En-Ma A 



(22) 



where the kinematical constants B = B' J [3(m A + m N )] and B' = ^2/3(M A /M N + {En(p) +M n ) /E N (p) and 
with kinematics where the A is produced at rest. Using Si and the other sink types used in our study of the electro- 
magnetic form factors l38l l42ll the four axial form factors C3 , C4 , Cj and Cg can be determined by performing an 
overconstrained analysis as described in Section IV. 

In Fig.0we show our lattice results for the four axial form factors for Wilson fermions and in the hybrid scheme. 
For the Wilson fermions we use Z A — 0.8 fislEfll . In all cases we observe that C3 is consistent with zero. Comparison 



XN r =0 m^O.56 GeV 

ON p =0 m^O.49 GeV 

*N p =0 rn^O.41 GeV 

AN t =2 m n =0.69 GeV 

• N F =2 m^O.51 GeV 

ON =2 m =0.38 GeV 




X N p =0 m^O.56 GeV 

O N,=0 m,=0.49 GeV 

* N,=0 m,=0.41 GeV 

* MILC m^O.59 GeV 

* MILC m^O.50 GeV 
■ MILC m n =0.34 GeV 



0.5 1.0 1.5 

Q 2 (GeV 2 ) 



0.5 1.0 

Q 2 (GeV 2 ) 



Q 2 (GeV 2 ) 



0.5 1.0 1.5 

Q 2 (GeV 2 ) 



FIGURE 7. The axial form factors C^C^, C% and c£ as a function of Q 2 . The notation is the same as that in Fig. [5] 
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of quenched and unquenched results obtained with Wilson fermions shows that unquenching effects are small for the 
dominant form factors, C5 and Cg . The form factor C4 shows an interesting behavior: The unquenched results for 
both dynamical Wilson and domain wall fermions show an increase at low momentum transfers. Such large deviations 
between quenched and full QCD results for these relatively heavy quark masses are unusual making this an interesting 
quantity to study effects of unquenching. 



For finite mass pions partial conservation of axial current (d^A^ix) — f n m 2 n a (x)) leads to the off-diagonal 
Goldberger-Treiman relation C^(Q 2 ) = fxgnN\{Q 2 ) /2Mn where gnNAiQ 2 ) is determined from the matrix element 

of the pseudoscalar density < A + | \\i{x)y^ Wi x )\P > an d me pi° n decay constant f n from the two-point function 
< 1 A4 (x) I K > . To relate the lattice pion matrix element to its physical value we need the pseudoscalar renormalization 
constant, Z p . For quenched lE5ll and dynamical Wilson fermions |46] Z p (fj, 2 a 2 ~ 1) = 0.5, which however may depend 
on the renormalization scale. At the chiral limit, at the lowest available momentum, we find g^NA 

(g 2 = 0.135GeV 2 ) = 

18.0(1 .9) and g K Nk{Q 2 = 0.443 GeV 2 ) = 15.8(1 .8) for quenched and dynamical Wilson fermions respectively where 
the quoted errors include a 10% uncertainty in Z p . These values are to be compared with gxNA 

(m 2 ) = 23.2 ±2.6 

obtained from an analysis of nN scattering |47]. We show in Fig. [H]the ratio fngnNA./ (2MnC$) for Wilson fermions. 
We observe that this ratio is almost independent of Q 2 and and that, as the quark mass decreases, it becomes consistent 
with unity in agreement with the off-diagonal Goldberger-Treiman relation. 

A prediction of lattice QCD is the ratio C5 jC\. The form factor C\ can be obtained from the electromagnetic N to 
A transition. Using our lattice results for the dipole and electric quadrupole Sachs factors, ^1 an d 5^£2> we extract C3 
using the relation 

v 3 M A (M N +M A ) 
C i=2(M N +M & ) 2 + Q 2 ^-^E2). (23) 

The ratio C5 /CY, shown in Fig. [S] for m n ~ 500 MeV, has values that fall within errorbars in the quenched theory 
and in full QCD. Given this agreement between quenched and unquenching results in this quark mass range, we opt to 
extrapolate the quenched results, which carry the smallest errors, to the chiral limit. As can be seen from Fig.[8]a linear 
extrapolation in ;« 2 leads to only a small increase in this ratio at the chiral limit. Without an analysis within a chiral 
effective theory for the quark mass dependence, this is the best that can be done to extract a reasonable estimate for the 
physical value of this ratio. Under certain assumptions, such as taking C3 ~ and considering that C4 is suppressed 
as compared to C5 , both of which are justified by the lattice results, the parity violating asymmetry can be shown to be 
proportional to this ratio |48]. Our lattice results show that this ratio and, to a first approximation the parity violating 
asymmetry, is non-zero at Q 2 = and increases by a factor of 2-3 when Q 2 ^1.5 GeV 2 . 



VI. CONCLUSIONS 

State-of-the-art lattice QCD calculations yield accurate results on a number of observables that are important for 
understanding the structure of the nucleon. We have presented the framework of such lattice computations for the 
nucleon and N to A transition form factors for quenched and dynamical Wilson fermions. The results are accurate 
enough to allow a meaningful comparison to experiment. We also evaluated the N to A transition form factors within 
a hybrid scheme, that combines the best simulation of the QCD vacuum that is available up to now using staggered 
fermions, with domain wall fermions the have good chiral properties. Within this hybrid scheme, we obtain results 
which are consistent with those obtained with dynamical Wilson fermions, albeit with larger statistical errors. Although 
chiral fermions are more expensive, smaller quark masses can be reached without any conceptual difficulties. Work is 
in progress to calculate the A form factors and coupling constants within the same framework. The aim of this program 
is to evaluate a complete set of observables for the nucleon-A system at small enough quark masses so that one can 
use only lattice input to fix the parameters of chiral effective theories. Extrapolation to the physical regime can then be 
carried out avoiding uncontrolled approximations. Dynamical Wilson and MILC configurations at pion masses of 250 
MeV on large enough volumes are now becoming available and therefore, using the technology developed, we will be 
able to calculate these quantities closer to the chiral limit where we can reliably make contact with chiral perturbation 
expansions. 
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